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INTRODUCTION 

Microparticles can deposit on surfaces under the in- 
fluence of a variety of forces. In the absence of grav- 
ity, forces that would otherwise be negligible become 
dominant. Even in a gravitational field some of these 
forces can be much larger than the gravitational force. 
Among these are the phoretic forces that arise when a 
small particle exists in a nonuniform gas. In the case 
of the thcrmophoretic force, molecules colliding with 
the particle that come from a higher temperature region 
exchange more momentum and kinetic energy with the 
particle than do molecules coming from a lower tem- 
perature region. The resulting force can cause the par- 
ticle to deposit on a colder surface or move away from 
a higher temperature surface. Applications include 
particle deposition on semiconductor surfaces and on 
surfaces in combustion processes, containerless proc- 
essing, and the production of nanophasc materials, 
pigments and other fine particles. 

Momentum and energy transfer associated with gas 
molecule collisions are not well understood and cannot 
be predicted from first principles for “engineering” 
surfaces. Neither of the two limiting cases of perfect 
accommodation or specular reflection adequately de- 
scribe most collisions, so in the last century James 
Clerk Maxwell introduced the concept of accommoda- 
tion coefficients to describe collisional processes. Ex- 
periments are needed to determine accommodation 
coefficients. 

This research program involves theoretical and ex- 
perimental studies of phoretic forces that act on mi- 
crometer-size particles. Phoretic forces are strongly 
dependent on the mean free path and kinetic energy of 
the gas molecules, and for a given microparticle 
monatomic, diatomic and polyatomic gases behave 
very differently. Li and Davis [1,2] reported meas- 
urements using helium, nitrogen and carbon dioxide 
for a variety of microspheres including nonvolatile 
liquid droplets, polystyrene latex (PSL) spheres, glass 
microspheres and metallic nickel spheres. These sys- 
tems represent a wide range of thermal conductivities 
and other properties. Li and Davis compared their 
results with numerous theories in existence, and to 
make such comparisons it is necessary to know the 
thermal properties of the gas and microparticle. Fur- 
thermore, their work was limited to spherical particles. 
Many microparticles of interest have irregular shapes 
and unknown thermal properties, so it is necessary to 
measure their thermal properties to interpret and pre- 
dict the thermophoretic force. 

To perform thermophoretic force experiments we 
have built a special electrodynamic balance (EDB) 


equipped with heat exchangers mounted in a vacuum 
chamber. The chamber pressure can be varied to alter 
the mean free path of the gas molecules. The thermal 
and radiative properties of microparticles can also be 
measured by single particle experiments in which a 
charged microparticle is levitated in an EDB and illu- 
minated with a laser beam. The principles and appli- 
cations of the EDB have been surveyed by Davis [3,4]. 
Spjut et al [5] at MIT introduced the electrodynamic 
thermogravimetric analyzer (EDTGA) for the study of 
carbonaceous particles, and the MIT group performed 
a variety of investigations reviewed by Bar-Ziv et al. 
[6]. These include particle temperature measurements 
by infrared pyrometry and transient heating experi- 
ments [7]. Bar-Ziv and Sarofim [8] provided an addi- 
tional review of the literature associated with particle 
heating in an EDB, including the effects of microparti- 
cle properties on the electromagnetic heat source func- 
tion. 

Monazam et al. [9] and Monazam and Maloney 
[10] al the Morgantown Energy Technology Center 
(METC) adapted the techniques of Spjut and his co- 
workers to determine heat capacities, temperatures and 
absorptivities of single carbon particles by means of 
pulsed heating experiments. A levitated particle with a 
diameter of order 100 pm was illuminated from two 
sides using a CO 2 laser (10.6 pm wavelength) with a 
repetition rate of 100 Hz and a pulse duration of 3 ms. 
The particle temperature was determined from meas- 
urements of the radiant emissive power of the particle 
using Planck's distribution law for the monochromatic 
radiant emission intensity. 

The time-dependent temperature distribution in the 
microsphere was modeled by assuming that the heat 
flux associated with the electromagnetic radiation is 
uniform over the surface of the sphere. That is, the 
heat source was treated as a boundary condition rather 
than as an electromagnetic heat source within the 
sphere. Convective and radiative heat losses were 
included in the surface boundary condition, but for the 
conditions of the experiments gas phase conduction 
was the dominant heat loss. They solved the govern- 
ing unsteady state heat conduction equation numeri- 
cally, and performed parametric studies, comparing 
predicted surface temperatures with measured tem- 
peratures. The authors recognized that the use of a 
spatially uniform energy flux at the surface is an ap- 
proximation that may not be valid for other materials 
with physical and optical properties that differ from 
carbon. 

It was the objective of part of our recent studies to 
develop a more rigorous analysis of unsteady state 
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heating of a sphere by electromagnetic radiation, tak- 
ing into account a spatially-dependent surface energy 
flux based on Mie theory. The analysis of the METC 
group is recovered as a special case of a more general 
formulation. 

PULSED HEATING OF 
MICROSPHERES 

The governing conduction equation and auxiliary 
conditions used by Monazam and his coworkers are 
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ditions of the METC researchers, h r was less than 10% 
of h c , so the linearization should be accurate at later 
times. The heat loss at early times is overestimated by 
this approximation because h r is overestimated by us- 
ing the asymptotic mean surface temperature. Thus, it 
can be anticipated that the analytical solution devel- 
oped using approximation (4) will underpredict the 
temperature at small times but should be accurate after 
some transient period. 

The solution of the linearized problem is 
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in which p, Bi, x and x are dimensionless parameters 
and variables defined by 


P = aa r I inc /2 k,T 0 , Bi =(h c +h r )a/ K, , 
x= r/a, x = K,t/ a 2 p,C p i . 


( 7 ) 


The boundary condition at the surface of the sphere 
includes the time-dependent intensity, I inc ,(t), of the 
pulsed incident beam, the convective heat transfer co- 
efficient, h c , and the radiative heat loss. The absorp- 
tivity, otr is assumed to equal the emissivity, 8^, of the 
particle. Here, is the thermal conductivity of the 
microparticle, and <7 is the Stefan-Boltzmann constant. 
The incident intensity in the surface boundary condi- 
tion is divided by two because the laser beam was split 
to heat the particle from two sides. 

Monazam et ai solved this system of equations 
numerically, but for the conditions used in their ex- 
periments an approximate analytical solution can be 
obtained. This is accomplished by writing the linear 
approximation 

T 4 (a,t) -T 4 =h r [T(a,t) — T«,], (4) 

in which h r is a radiant energy heat transfer coefficient 
defined by 

h r =[Tm(a,-) + TijT m (a,oc)+T 00 ] (5) 

and T m (a,oo) is the mean surface temperature at large 
times, which can be computed by solving an energy 
balance on the sphere under the assumption that the 
heat flux averaged over one cycle is zero at large times 

The heat transfer coefficient, h c , is K 2 /a for a stag- 
nant gas and for low gas Peclet numbers. For the con- 


Here the subscript 1 refers to the particle and 2 denotes 
the gas, P is the dimensionless electromagnetic energy 
flux, Bi is a composite Biot number, and f(x) is the 
normalized laser pulse function, which is unity for the 
first 3 ms of a pulse and zero for the next 7 ms for the 
experimental conditions at METC. 

The eigenvalues in Eq. (5) satisfy the transcenden- 
tal equation 

tanX n = /(] -Bi). (8) 

The solution represented by Eq. (6) is compared 
with the numerical solution of Monazam and Maloney 
in Figure 1 . The physical parameters they used are a = 
70 jam, otr = Er = 0.85, line = 4.50 MW/m 2 , Kj = 1 .675 
W/m-K, k 2 = 0.048 W/m-K, C\ = 1466 J/kg-K, and pi 
= 1050 kg/m 3 . The Biot number corresponding to 
their conditions is 0.032. The agreement is excellent 
at later times, and during the transient period when the 
time-averaged temperature changes, the analytical so- 
lution predicts temperatures only slightly lower than 
the numerical solution. 

A more detailed comparison between theory and 
experiment is shown in Figure 2 for one cycle in the 
asymptotic region. The agreement is very good, indi- 
cating that the best-fit parameters selected by 
Monazam and Maloney are consistent with their model 
of uniform surface heating. 
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Figure 1. A comparison between theory and experi- 
ment for a carbonaceous particle. 



Figure 2. A comparison between theory and ex- 
periment for one cycle at later times. 


SPATIALLY NONUNIFORM HEATING 

The assumption of uniform surface heating is 
highly questionable when the thermal conductivity of 
the particle is not as large as that of carbon. For the 
140 pm diameter sphere of Monazam and Maloney, 
the light scattering size is X =2na/X= 41 .498. Using a 
refractive index of iV= 5.0+i4.0, we computed the 
source function, EE*/E 2 inc , presented in Figure 3 for 
one-sided heating using Mie theory as outlined by Bo- 
hren and Huffman [9]. Due to the strong absorption of 
electromagnetic energy, the source is concentrated 
near the surface, as assumed by Monazam and Ma- 
loney, but the source is highly nonuniform in the 0- 
direction. Our computations indicate that all of the 
energy absorption is confined to the outer 2% of the 
sphere, that is, to the region x = r/a > 0.98, and the 


source function is nearly independent of azimuthal 
angle <|). 


o 

tr 
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Figure 3. The internal heat source function for a car- 
bonaceous sphere with X = 41 .498 and N = 5.0+14.0. 



The internal heat generation function for an ab- 
sorbing sphere is 
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The pulsed heating problem can now be solved 
using a source function such as that depicted in Figure 
3 or, alternately, as an equivalent surface source prob- 
lem. For two-sided heating, the shape of the source 
function, considered as the surface source shown in 
Figure 3, can be approximated by the intensity distri- 
bution 


l(a,0)= I max cos 2 0. (10) 

Here 1,^ is the intensity at 0 = 0°,180°, and 0 is 
measured from the direction of propagation of one of 
the laser beams. This approximation is in very good 
agreement with the exact surface intensity distribution 
computed from Mie theory as demonstrated in Figure 
4. 

When the angular dependence of the temperature 
distribution is considered, the nondimcnsional heat 
conduction equation becomes 

dx x 2 dx 1^ X 3x ) x 2 3 t| 
in which 

U = ” — , and rj = cos0. (12) 
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Figure 4. Comparison between the exact solution for 
the surface intensity distribution and Eq. (10). 


The auxiliary conditions are given by 
U(x,T|,0)= (0,11, t) = 0, (13) 

OX 

and 
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The solution for the surface temperature is 
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Here N m0 and N m2 are the norms of the eigenfunctions 
X m0 (x) and X m2 (x), respectively, given by 

N m0 = t/(y m0 -sinYmQCOSYmo)/ JT/mO . (16) 


and 

N m2 = 

J— ^5/2(Ym2) + J3/2(Ym2)f“ 2 y m2 J 5/ 2(Ym2 ) J 3/ 2(Ym2 )■ 

(17) 


The eigenfunctions may be written in terms of Bes- 
sel functions as 


X m0 (x)=7;,Ji/2(Ym0x). O 8 ) 

and 

X m 2( x ) = J 5/2 (Y m2 x )• 0^) 

The eigenvalues y m o satisfy Eq. (8), and the eigen- 
values Y m2 satisfy 


tany m2 = 
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The integrals I u I 2 and I 3 are defined by 


I, =Jx 2 X m0 (x)dx, 

(21) 
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For Biot numbers larger than those for relatively 
highly conducting carbon spheres Eq. (15) predicts 
very large angular variation in the surface temperature. 
Such variations must be taken into account in the in- 
terpretation of data to determine the thermal properties. 
The surface temperatures at 0 = 0° and 90° arc illus- 
trated in Figure 5 for K } =0.167 W/m-K (Bi = 0.320). 
The other physical properties are those used by 
Monazam and Maloney. Additional results have been 
reported in a paper by Davis and Widmann [12] that is 
in press. 
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Figure 5. Surface temperature histories for a lower 
thermal conductivity particle. 
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THERMOPHORET1C FORCE 

The existing apparatus used for thermophoretic 
force measurements, which is shown in Figure 6, has 
been modified to improve the data acquistion system 
and to facilitate the alignment of the electrodes used 
for particle levitation. Studies of the effects of natural 
convection on the thermophoretic force have been car- 
ried out over a wide range of system pressures, and 
Figure 7 shows the levitation voltage as a function of 
pressure for a polystyrene latex sphere in air. 

The change in the levitation voltage (with and 
without a temperature gradient in the gas phase) yields 
the ratio of the thermophoretic force to the gravita- 
tional force. If V 0 is the dc levitation voltage required 
to balance the gravitational force, mg, and V is the 
levitation voltage required to balance the particle when 
gravity and other external forces, F, act on the particle, 
the force ratio F/mg is given by 
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Figure 7. Levitation voltage data for a PSL sphere in 
air as a function of the chamber pressure. 


F / mg = (Vo - V)/ V 0 . (24) 



Figure 6. Cross section of the electrodynamic balance 
used for thermophoretic force measurement. 


The external forces can include the phoretic force 
as well as aerodynamic drag associated with convec- 
tive motion of the gas. 

The data of Figure 7 indicate that for pressures 
greater than 20 torr natural convection occurred due to 
edge effects associated with temperature differences 
between the heat exchangers and the chamber wall, but 
those effects become negligible as the system pressure 
decreases. At low pressures the free-molecule regime 
is reached, and the force becomes independent of pres- 
sure. 


Data for the intermediate or transition regime (be- 
tween the continuum regime and the free-molecule 
regime) are used to determine the accommodation co- 
efficients for momentum and energy transfer between 
the gas molecules and the microparticle.lt is highly 
desirable to minimize the effects of natural convection, 
and modifications in the apparatus will be made to do 
so. 
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